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Abstract 

Solitons in the Skyrme-Faddeev model on M 2 x S* 1 are shown to 
undergo buckling transitions as the circumference of the S l is varied. 
These results support a recent conjecture that solitons in this field 
theory are well-described by a much simpler model of elastic rods. 

1 Introduction 

It has recently been conjectured [Ij that solitons in a particular field 
theory, the Skyrme-Faddeev model [2], are well-described by an effec- 
tive model based on Kirchhoff elastic rods. It was shown in [T] that 
the elastic rod model gives a good qualitative approximation and a 
reasonable quantitative approximation to low-charge Skyrme-Faddeev 
solitons. 

The discovery of this elastic rod model motivates the search for 
elastic phenomena, such as buckling, in the Skyrme-Faddeev model. 
In the present article we will investigate one such buckling effect in 
both the elastic rod and the Skyrme-Faddeev models, caused by the 
simultaneous stretching and twisting of a length of elastic rod (or 
soliton). On a technical level, the most convenient way to stretch 
elastic rods and Hopf solitons is to place them on the manifold IR 2 x S 1 . 
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The rod is arranged to wind once around the S , and can be stretched 
by varying the circumference P of S . 

The kinds of effects that can occur are sketched in figure [T] The 
simplest configuration consists of a straight rod (a). If it is twisted in 
a suitable way, this straight rod may buckle at some critical value of 
P to form a helix (b). This helix could buckle again at another value 
of P, forming a kinked configuration (c). 




Figure 1: Sketch of helical buckling 

The straight rod (a) is fixed by the group SO(2) x SO(2), where 
one copy of SO (2) acts on S 1 and the other by rotation on M 2 . The 
helix (b) is invariant only under a diagonal subgroup SO(2), while the 
kinked configuration has no continuous symmetries. Thus the buckling 
transitions are examples of spontaneous symmetry breaking. 

The motivation behind our research is twofold. On the one hand, 
we were interested to further test the reliability and utility of the 
elastic rod model as a description of Skyrme-Faddeev solitons. The 
rod model has already been successfully used to describe solitons on 
M 3 , and R 2 x S 1 seems a good place to test it further. 

On the other hand, we wanted to explore what seem to be fairly 
important aspects of the Skyrme-Faddeev model. Solitons in this 
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model resemble knotted rods, so it seems a good idea to investigate 
basic properties of these rods, such as their response to stretching and 
twisting. Some studies of the Skyrme-Faddeev model on R 2 x S 1 have 
appeared before [3l 0], but, surprisingly, the simple buckling effects 
described here have not previously been investigated. The Skyrme- 
Faddeev model has been proposed as a model of glueballs [5]. 

An outline of the rest of this article is as follows. In section [2] we 
will review the Skyrme-Faddeev and elastic rod models on R 3 , and the 
connection between them. In section [3] we explain our conventions for 
putting these models on R 2 x S , and discuss in detail the topological 
charges of Skyrme-Faddeev solitons on this space. In section [4] we 
present some analytical and numerical methods that can be used to 
study the elastic rod model on M 2 x S 1 . In section [5] we present our 
main results, including direct comparisons between the two models. 
We draw our conclusions in section [6l 

2 The Skyrme-Faddeev model and elas- 
tic rods 

2.1 The Skyrme-Faddeev model 

The Skyrme-Faddeev model is an 0(3) sigma model in 3+1 dimen- 
sions, whose Lagrangian is augmented by an additional term quar- 
tic in derivatives. We will only be interested in static states, so 
that the fields of the model can be taken to be a triplet of scalars 
<ft a = (cj) , cj) 2 , cj) S ), which are functions of x G M 3 , and which satisfy the 
constraint <p a (p a = 1. The static theory is defined by specifying the 
energy functional, 

where 

Fij = e^rdi^d^ . (2) 

Static states in the model are solutions of the Euler-Lagrange equa- 
tions associated to E$f- 

A configuration <j> a {x) has finite energy only if di4> a — > as |x| — > 
00. Hence, for finite-energy configurations the limit of cp a (x) as |x| — > 
00 is a well-defined point on the 2-sphere. Therefore finite-energy con- 
figurations can be extended to continuous maps from S 3 = M 3 U {00} 
to S 2 . It is well known that such maps have a topological invariant 
Q € 7Ts(S 2 ) = Z, known as the Hopf degree or topological charge. 
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The Hopf degree can be calculated in one of two ways. The first 
possibility involves looking at the preimages of points on S 2 under the 
map cf) a . Generically these preimages will be unions of disjoint loops in 
IR 3 . The Hopf degree Q is equal to the linking number of the preimages 
of two distinct points [6]. Alternatively, Q may be calculated using 
an integral formula. For any finite-energy configuration the tensor F^ 
defines a 2-form on the 3-sphere. This is closed, and also exact since 
H 2 (S 3 ) = 0. Therefore one can find a 1-form A such that F = dA, 
and Q is equal to the integral, 

Q = vhl AAK (3) 

The most important problem in the Skyrme-Faddeev model is the 
identification of stable static states. Numerical simulations indicate 
that for each value of Q there exists a unique energy-minimising con- 
figuration. It is known that the energies of these configurations scale 
like Q 3 / 4 ; more precisely, it was proven in [71 [8] that there exist con- 
stants C\ and C*2 such that 

CiQ 3/4 < inf E SF < C 2 Q 3/ \ (4) 
Q 

where the infimum is taken over fields with Hopf degree Q. Conjec- 
turally, the value of the constant C\ can be taken to be 1 [9]. 



2.2 Elastic rods 

In [lj it was shown that stable static states in the Skyrme-Faddeev 
model are well-approximated by elastic rods. An elastic rod consists 
of two vector- valued functions x, m of a real parameter a. The first 
function x specifies the location of the centreline of the rod in M 3 . The 
second function fn specifies how the rod is twisted, and must satisfy 
the constraints, 

m(a) ■ fh(a) = 1, fh(a) ■ x'(a) = 0. (5) 

For small e, the function x(a) + em(a) describes the location in M 3 of 
a straight line in the material of the rod close to the centreline. 
The energy functional for elastic rods is En = AL + Ek , where 

L = J \x'\da (6) 

is the length of the rod and 

E K = [ (Bk 2 + Cj 2 ) \x'\da (7) 
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is Kirchhoff's energy functional. Here k is the curvature of the rod, 
defined by 

^) = iS' (8) 



where 



is the unit tangent vector to the rod. On the other hand, 7 is the 
twist rate of the rod, and is defined by 

, . t ■ mf x m 1 . 

l{a) = . (10) 

By construction, the elastic rod energy functional is independent of 
the parametrisation. It is often convenient to choose a parametrisation 
for which \x'\ = 1, in which case the parameter a is denoted s and 
called the arclength parameter. The energy functional is invariant 
under SO(2) rotations of the normal vector m, generated by 

5fh = t x m. (11) 

A slightly different way of describing elastic rods involves the Frenet 
frame. The Frenet frame can be defined when the centreline x(s) is 
arclength-parametrised, with arclength parameter s. It consists of 
three vectors t(s), n(s) and 6(s), where t is the unit tangent vector, 
n = t /k, and b = t x n. These vectors satisfy the Serret-Frenet 
equation: 

A ( t \ ( k \ / t \ 

(12) 

where r(s) is a real function known as the torsion. As long as k / 
the material frame vector of a rod can be written m = sin an + cos a b 
for some real function a(s). Then the twist rate is 

7 (s) = a'(s)-r(s). (13) 

In order to obtain a better match with Hopf solitons, we will impose 
a non- intersection constraint on our elastic rods. We will assume that 
the rods have a circular cross-section of radius p, and demand that the 
rods do not intersect themselves. According to [10], this is equivalent 
to the following two conditions: 

1. k(s) < p^ 1 for all s, 
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2. I > p, where 

I = -min {d{s\, S2)\s\ / S2 and (s\, S2) is a critical point of d} 

and d(si,S2) is the distance between the points x(s\), x(s2)- 

The first of these conditions essentially says that the radius of curva- 
ture k _1 of the centreline cannot be less than p. 

2.3 Elastic rods from Skyrme-Faddeev solitons 

The purpose of this article is to compare energy minima in the Skyrme- 
Faddeev and elastic rod models. In order to make this comparison, we 
need to explain how configurations in the two models are related. We 
do this by specifying a map from finite-energy field configurations in 
the Skyrme-Faddeev model to configurations of elastic rods. A map 
sending elastic rod configurations to Skyrme-Faddeev field configura- 
tions was described in pp. 

First of all, the centreline of the rod can be defined to be the 
preimage under (j) a of a point in S 2 antipodal to the asymptotic value 
of 4> a . More concretely, by making an SO(3) rotation of <p a {x) we can 
arrange that the following boundary condition is satisfied: 

(j) a (x) -> (0,0, 1) as r -»• 00. (14) 

Then the centreline of the rod is defined to be the collection of all 
points x such that <j) a (x) = (0,0,-1). Typically this set will consist 
of a number of closed loops, each of which can be parametrised as 
xj(<ti) with / an index labelling the loop. Since the loops are closed, 
the parameter 07 can be chosen to lie in a closed interval [<7j,0j], in 
such a way that 

x I (a° I ) = x I (a}). (15) 

The material frame vector fhj(aj) is obtained by projecting the 
vector 

(16) 

orthogonally onto the space perpendicular to x'j(ai), and normalising. 
The obtained functions rhj(ai) satisfy 

m 7 K) = rnjia}). (17) 

Notice that the SO (2) rotations of the material frame vector m corre- 
spond to the SO (2) C SO (3) which fixes the asymptotic value of 4> a . 
We define the charge Q of a collection of elastic rods to be the link- 
ing number of the collection of curves xi with the collection of curves 
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xj + em j, for small enough e. Then the map from Skyrme-Faddev 
field configurations to elastic rod configurations obviously preserves 
Q. 

In [1] it was shown that, for suitable choices of the parameters 
A,B,C,p, minimum-energy configurations in the elastic rod model 
with some fixed value of Q look similar to minimum-energy configu- 
rations in the Skyrme-Faddeev model with the same value of Q. The 
match was obtained by choosing the dimensionless parameter C/B to 
be 

C/B = 0.85, (18) 
and fixing the rod thickness p to be 



The remaining two parameters correspond to choices of units of length 
of energy. We will choose these so that the energy minima in the 
elastic rod model with 1 < Q < 7 have similar sizes and energies to 
the corresponding solitons in the Skyrme-Faddeev model, as in [T]. 
This leads to 

A = 0.0872, B = 0.0671, C = 0.0571. (20) 

These parameters give the charge 1 soliton the correct size. They over- 
estimate its energy (by about 10%), and underestimate the energies 
of solitons with higher charge. 



3 Solitons wrapping a circle 

We will study Skyrme-Faddeev solitons not on R 3 , but on the space 
K 2 x S 1 . The configuration space is by definition the set of <S 2 -valued 
functions 4> a (x) satisfying 

<f> a {x l ,x 2 ,x 3 + P) = (f> a (x\x 2 ,x 3 ), (21) 

for some P > 0. Static stable states will be local minima of the energy 
per period, 

E SF = / (d^r + \F tJ F^)d 3 x. (22) 

Configurations with finite energy per period must satisfy the boundary 
conditions, 

d r 4> a -> as r -»• oo, (23) 
d z 4> a -> as r -> oo, (24) 
d d (j) a -> as r ^ oo, (25) 
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where r, 9 are polar coordinates on M. and z = x is the coordinate in 
the periodic direction. These boundary conditions imply that <j) a has 
not one but two conserved topological charges, as we now explain. 



We begin by considering the boundary condition (23). This implies 
that 4> a has a well-defined limit as r 
may depend on 9 and z. 



oo. The limiting function 



a 

oo 



The boundary condition (24) implies that <^ does not depend 
on z, but may still depend on 9. It follows that the map 4> a can be 
extended to the compactification of M 2 x 5 1 obtained by adding a 
circle parametrised by 9 at infinity. This compactified space is in fact 
the 3-sphere. One way to see this is to consider the following map 
from R 2 x S 1 to S" 3 : 



/ 



(r, 9, z 



\ 



(26) 



The reader may verify that this is an injection, and that the boundary 
of the image of M. 2 x S 1 is a circle parametrised by 9. So, any field 
configuration 4> a satisfying ( 24 ) extends to a map from S* 3 to S 2 , and 



hence has a Hopf degree QgZ. The Hopf degree may be calculated as 
above, either by determining the intersection number of the preimages 
of two distinct points, or by integrating A A F over R 2 x 5 1 . 



Now consider the boundary condition (25). This condition implies 
that cj)^ does not depend on 9, but does not rule out z-dependence. 
It follows that <f> a can be extended to the compactification of M 2 x 5 1 
obtained by adding a z-dependent circle. This compactified space is 
S 2 x S , since the 1-point compactification of M 2 is S 2 . Maps from 
S 2 x S 1 to S 2 have two topological charges [HJd!]. First, there is a 
charge DgZ, which is equal to the degree of a map from S 2 to S 2 
obtained by restricting to a slice of constant z (this is independent of 
the choice of slice). And second, there is a Hopf degree Q' eZjj. 



Finite-energy configurations satisfy all three conditions (23), (24) 



and (25), so have three charges Q, Q' , D. These three charges are not 



independent: Q is equal to the value of Q modulo D. Therefore there 
are two independent topological charges, Q,D G Z. 

In the present article we will restrict attention to configurations 
with D = 1. The corresponding elastic rods will have only 1 strand 
and will consist of functions x, rh of a G M, satisfying 



x{a + a P ) = £(a) + (0,0,P), fn(a + a P ) = m(a), (27) 
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for some op > 0. The energy per unit period is 

rap 

E R = (A + Bk 2 + Cj 2 )\x\da. (28) 
Jo 

Here we will investigate Q = 1,2,3; Skyrme-Faddeev solitons with 
D = 1 and 5 < Q < 8 have previously been studied in f3j. Skyrme- 
Faddeev solitons with D > 1 have been investigated in and these 
correspond to elastic rods with D strands. 



4 Buckling of elastic rods 

In this section we will outline numerical and analytical methods that 
can be used to study the buckling of elastic rods on M 2 x S 1 . We begin 



by considering in subsection 4.1 a straight rod winding around the S . 
In subsection 4.2 we will consider a more general helical configuration, 
and derive conditions which determine when the straight rod may 
buckle to form helix. In subsection|4.3|we briefly consider how the non- 



intersection constraints may be applied to the helix. In subsection 4.4 



we study buckling of the helix itself, and in subsection 4J5 we describe 
some numerical methods for studying elastic rods. 

4.1 Straight rod 

We begin by making an ansatz: 

x(s) = (0,0,s) , . 

m(s) = (cos(2ttsQ/P), - sin(27rsQ/P), 0). 



This ansatz satisfies the boundary conditions (27) for elastic rods. 
The Hopf degree is Q and the second topological charge is D = 1. 
The ansatz has an SO(2) x SO(2) symmetry. The first copy of SO(2) 
acts by translation on z and rotates rh, and the second copy acts by 
rotation on x, y and also on fh. In fact, the ansatz is the unique ansatz 
with these symmetries, so by the principle of symmetric criticality it 
solves the equations of motion for elastic rods. The energy per period 



of the ansatz ( 29 ) is 



E R = AP + ^. (30) 
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4.2 First buckling 

Now we consider a more general ansatz, which describes a helix with 
M coils per period P: 

/ . 2-kMs „ 2vrMs Ps 
x{s) = Rsm — - — ,/tcos 



L ' L ' L 

m(s) = sina(s) n(s) + cos a(s) 6(s) (31) 
2vr(Q - M)s 



as 



Here R > is a parameter describing the radius of the coils of the 
helix, and L is the length of the rod, given by the formula 

L 2 = P 2 + (2-kMR) 2 . (32) 

The reader may verify that \x'\ = 1, and hence that the rod is arclength- 
parametrised. 

The Hopf degree of the configuration (31 ) is independent of R, and 
is most easily evaluated in the limit R — )■ 0. First of all, the vectors 
n, b can be calculated to be 



/ 2ttMs 2vrMs 
n(s) = I — sin — — — , — cos — — — , 

. (P 2ttMs P 2-kMs 2ttMR\ 

b{s) = {l cos ^T>-l sm ~L-' — Z ) 

In the limit R — > these expressions reduce to 

/ 2vrMs 2vrMs 
n(s) = I -sin - , -cos - ,0 

/ 2vrMs 2ttMs 
b(s) = I cos— — , -sm— — ,0 



(33) 



(34) 



Therefore in the limit R — > the ansatz (31) reduces to 
x(s) = (0,0, s) 

-* / \ • -* i r ( 2vrsQ . 2vrsg \ (35) 

m(sj = sm a n + cos a o = I cos — — — , — sm — — — , 



which coincides with the ansatz (29) for a straight rod. In this limit 
the Hopf degree is obviously Q. 

For R > 0, the symmetry group of the ansatz (31) is an SO(2) 
subgroup of SO (2) x SO (2). When R = @) reduces to ((29]). So 
the helix can be regarded as a symmetry-breaking perturbation of the 
straight rod. 
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It is convenient to introduce a dimensionless parameter A = Pj L 6 
(0, 1]. The radius and length can be recovered from A using the for- 
mulae L = P/X, 2nMR = PVA 3 ^ 
helix is 



1. The energy per period of the 



+ A-1 



(36) 

When A = 1 this reduces to the energy (30) of the straight rod. We 
know that the straight rod is a critical point of the energy functional, 
but is it stable to small helical perturbations? 

To answer this question, we just need to look at the slope of the 
function Er{\) near A = 1. We have 



dE R 



dA 



ir 8ir 2 M 2 B 4vr 2 C^ 2 . 
-AP + -^{Q 2 + 2QM). 



A=l 



P 



P 



(37) 



The straight rod is unstable to the perturbation specified by M if and 
only if this derivative is positive, or equivalently, 



A (£) < C{Q 2 + 2QM) -2BM 2 . 



(38) 



The right hand side of this inequality does not depend on P. If the 
right hand side is positive one can always find non-zero values of P 
which satisfy the inequality, but if the right hand side is negative or 
zero the inequality can never be satisfied. Therefore this inequality 
can be satisfied if and only if the right hand side is positive. The right 
hand side is positive if and only if 




< M < 




(39) 



This means that for fixed Q only finitely many values of M give rise 
to instabilities of the straight rod. For example, when C/B = 0.85 
and Q = 1, 2, 3, only the following values of M need to be considered: 



(40) 



The value P c of P at which the straight rod becomes unstable to a 
helical perturbation is given by 



Q 


M 


l 


1 


2 


1,2 


3 


-1,1,2,3 



2?r 



CQ 2 + 2CQM - 2BM 2 



A 



(41) 
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0, therefore if the straight rod 



The energy (|36|) blows up as A 
at A 

some value Aq £ (0, 1). Since the ansatz (|3 1 h is fixed by symmetries 



1 is unstable, the energy (30) must attain a minimum at 
(0,1) 



this minimum of Er(X) corresponds to a solution of the equations of 
motion for the elastic rod model. 



4.3 Non-intersection constraint 



As discussed in section 2.2 demanding that an elastic rod does not self- 
intersect imposes two conditions on the configuration space of elastic 



rods. For helical rods (31), the first condition k < p 1 is equivalent to 



A 2 {1 - A 2 ) < 



(—)' 
\2irMpJ 



(42) 



The left hand side of this inequality is bounded above by 1/4, so if 
P/M > up this constraint does not restrict the range of A. If on 
the other hand P/M < np the constraint means that the allowed 
range of A is divided into two disjoint pieces (0, A_] U [A + , 1], where 
< A_ < A+ < 1 are defined by 2A| = 1 ± y/l - (P/Mvp) 2 . 

To understand the implications of the second constraint I > p, we 
first evaluate the distance function: 



d 2 (0,s) 



2-kMs 2ttMs Ps 

R cos — - — , it sin — - — , — 
Li Li Lj 



;i,o,o) 



p 



D(a) 



1 



2ir 2 V A 2 



1 



1 (1 - cos27rcr) +cj 2 . 



(43) 
(44) 
(45) 



The function D(a) may or may not have a local minimum a c S [0, 1], 
depending on the value of A. If a local minimum a c exists, the second 
constraint is satisfied if and only if 



D(a e ] 



(46) 



If on the other hand D has no local minimum, the constraint is sat- 
isfied. For values of A close to 1, D does not have a critical point, so 
a helix with A close to 1 always satisfies the second constraint. For 
values of A close to 0, D does have a local minimum a c . 

When it exists, the value of D(a c ) is less than 1, and tends to 1 
as A — > 0. This means that if P/M < 2p, helices with A ~ are 
completely ruled out (although helices with Awl are still allowed). 
This result matches geometric intuition: if P/M is less than twice 
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the thickness of the rod, a helix with large radius R cannot avoid the 
overlapping of neighbouring coils. 

In practise, the first constraint is more important than the second. 
For large values of P/M the energy function favours configurations 
with A ~ 1, so neither the first nor the second constraint influences 
the shape of the energy- minimiser. If P/M < np the first constraint 
may influence the shape of the energy-minimiser: in particular, there 
will be two local energy minima, since the first constraint divides the 
range of A into two pieces. The second constraint begins to affect the 
shape of the rod at even smaller values of P/M. 



4.4 Second buckling 

We have described above how the straight rod, with SO(2) x SO(2) 
symmetry, can buckle to form a helix, with only SO (2) symmetry. 
In this subsection we will describe how the helix can buckle again to 
form a kinked configuration with completely broken symmetry. We 
will present two tools with which this second buckling can be studied: 
first, we will present an analytical method for calculating the critical 
period P c at which the buckling occurs; and second, we will describe 
some numerical methods with which the kinked configuration can be 
studied. 

The analytical approach to the buckling is based on the treatment 
of Mitchell's instability of circular rods presented in |13j . It can be 
shown that the Euler-Lagrange equations for the elastic rod energy 
functional are 



= - A + Bk 2 + C7 2 + 2B[ r 2 I - 2C 7 r 



=B(2k't + kt') + Cjk' 
=7. 



(47) 



Here a prime denotes differentiation with respect to the arclength 



parameter s. The helix (31), with energy Er(\) given in (36), is a 
solution if and only if A solves the equation dE^/dX = 0. We will 
assume that A has been so chosen, and will ignore the restrictions 
imposed by the non-intersection constraints. The curvature, torsion 
and twist rate for the helix are 

27rM r 

k =— — vi - y 

r = - —A (48) 
2tt 

7o= T (Q + M(A-l)). 
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Now we suppose that a small perturbation of the helix has been 
made, so that k = kq + 6k, r = to + St, 7 = 70 + 6j. We will assume 
for simplicity that 



L i-L 

6k ds = / 6t ds = 0. 
Jo 

The linearised equations of motion are equivalent to 
°-fdJ-o + (^ + 2,o)>., 



(49) 



(50) 



with 6t = (Cjo/B — 2t)6k/k and 6j = 0. The linearised equations 
of motion have a solution 6k oc sin(27rns/L) only if there exists an 
integer n such that 



n 



M 2 (l 



A" ) + ( ^(Q + MX-M) 



2MXJ 



(51) 



The existence of a solution to the linearised equations of motion indi- 
cates the presence of a buckling instability. Thus it is straightforward 
to determine when buckling occurs: for each value of P one computes 



the value X(P) of A which minimises the helix energy (36), and from 



this, the right hand side of (51). Buckling can occur at any value P c 



of P for which the right hand side of (51 ) is the square of an integer. 



The above discussion was based on the assumption (49). It can 



be shown with a little more work that dropping this assumption does 
not yield any additional solutions to the linearised equations of mo- 
tion, essentially because perturbations for which 6k, 6t are constant 



correspond to modifications of the parameters P, X in the ansatz (31 ). 



4.5 Numerical methods 

After the helix has buckled, it is no longer possible to obtain analytic 
solutions for elastic rods. Instead, we employ numerical methods. A 
discretisation of the elastic rods was presented in [14J, which we briefly 
recall here. 

The centreline of the rod x(a) is replaced by a sequence of points 
x 1 with i E Zjy. The vector connecting two points is := x J+i — x 3 , 
and the length of this vector is lj = | e*,- 1 . The material frame m(cr) is 
replaced by a sequence of vectors fhj of unit length satisfying rhj ■ ej = 
0. 

A unit tangent vector is defined by tj := ej/lj. We also introduce 

Sl : = 1 (52) 



14 



The vectors fi J approximate n{(j)b{<j). The discretised curvature is 
calculated from 

{K j f = \W\ 2 . (53) 

Finally, a discretisation of the twist rate j(s) is given by 7*, defined 
via 

sin 



•fhj-i -tj-i x fl?)((V ■ rhj) — (rhj ■ tj x £V)($F ■ rhj-i) 



u j r 2 



(54) 



The discretisation of the elastic rod energy functional is now obvi- 
ously 

N ~ l 1+1 
E R =Y. A h + (B{^? + C(^) 2 )^^. (55) 

j=0 

It can be shown, with some effort, that this reduces to the usual energy 
functional in the continuum limit. We searched for minima of this 
discretised energy using a simulated annealing algorithm. The non- 
intersection constraint was imposed using the obvious discretisations 
of conditions 1 and 2 described in section 2.2 It was necessary to 
enforce the arclength parametrisation condition ^ ~ L\/i,j, in order 
to maintain a good approximation throughout the simulation. This 
was achieved by adding a penalty function to the energy. 



5 Skyrme-Faddeev solitons 

In this section we will apply the methods of the previous section to 
study in detail elastic rods on R 2 x S 1 with Hopf degree Q = 1, 2, 3. 
We then compare these predictions with full numerical simulations of 
the Skyrme-Faddeev model. 

Our simulations of the Skyrme-Faddeev model were performed by 
evolving the models equation of motion on a discrete lattice of 300 x 
300 x P/0.08 points, with lattice spacing Ax = 0.08. By varying the 
number of lattice points and the lattice spacing this size lattice was 
found to give the lowest energy solutions for all charges and periods P. 
The spacial derivatives were calculated on the lattice to fourth order. 
A Lagrangian multiplier was also included to preserve the constraint 
cj) a (j) a = 1. 
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5.1 A straight Skyrme-Faddeev soliton 

Before discussing full numerical simulations, we first consider a straight 



Hopf soliton analogous to the straight rod ( 29 ) . Consider the following 
ansatz: 



■ tt \ (a^ 2 ^ z \ ■ ft \ ■ (o^Q z \ ft \\ 
sin / (r) cos I U H — — I , sm / (rj sm I v -\ — — I , cos j(r) I . 

(56) 

This ansatz is invariant under a certain action SO(2) x SO(2), and is 
the most general ansatz with this symmetry. Substituting this ansatz 



into the energy functional (22) gives 



£ — I^f( ( « 2 + (? + ^)( 1 + ( « 2 ) si " 2 ^ dr - 

(57) 

We impose the boundary conditions, 

/(0) = vr, f(r) -> as r oo, (58) 

so that the ansatz is well-defined at r = and can have finite energy 
per period. Then the topological charge per unit period is Q. 



For any values of Q,P, the minimum of (57) with respect to vari- 
ations in / can be determined numerically, using a gradient flow al- 
gorithm. In fact, it is sufficient to do this for Q = 1, since the energy 



density in (57) only depends on the ratio Q/P. 



5.2 Charge 1 

In figure [2j we have plotted the energies of elastic rods with Hopf 
degree Q = 1 as a function of the period P. For large enough periods, 



the only state in the elastic rod model is a straight rod (29). As 
the period P decreases, the energy of this state decreases, attains a 
minimum at P = 5.08, and begins to increase. At P = 4.09 the 



straight rod becomes unstable and buckles to form a helix (31) with 
M = 1. 

At P = 3.75, just after the helix has formed, the non-intersection 
constraint begins to play a role. The space of allowable A is split 
into two intervals, and accordingly the elastic rod energy has two 
local minima, one being a helix with large radius R and the other 
being a helix with small radius. For most periods the helix with small 
R has the lowest energy. For a small range of periods [2.59,3.19] 
the helix with large R has the lowest energy. However, when P = 
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N = 1 



soliton 
straight soliton 
straight rod 
M=1 helical rod (small R) 
M=1 helical rod (large R) 



period 



10 



Figure 2: Energies of solitons and rods with Q = 1 

2.59 is reached the helix with large R is no longer permitted by the 
non-intersection constraint, and once again the small-i? helix is the 
favoured configuration. 

The main qualitative predictions of the elastic rod model are that 
for large periods the favoured configurations is a straight soliton, and 
for small periods the favoured configuration is a soliton in the shape 
of a helix with small R. These states have been depicted in figure 
[3j with a yellow tube representing the curve x(s) and a red tube 
representing the curve x(s) + em(s) for small e. There may also exist 
a helical soliton with large R for a small range of periods, but since 
the corresponding state is short-lived in the rod model, one cannot be 
confident that it would exist in the Skyrme-Faddeev model. 

In figure [4] we have shown pictures of the Skyrme-Faddeev solitons 
with Q = 1. In these pictures, the blue surface represents the preimage 
under the map <ft a : M 3 — > S 2 of a circle surrounding the south pole 
(ft 3 = —1 in S 2 , and the red surface represents the preimage of a circle 
surrounding a point near the south pole. The red tube links once with 
the blue tube, confirming that the charge is 1. The solitons appear at 
first sight to be straight, however for smaller periods a slight buckling 
can be detected. The pictures look similar to to the rods depicted in 
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Figure 4: Solitons with Q = 1 and P = 7, 5, 3, 1. 



figure |3j The main difference occurs at P = 3, where the buckling of 
the rod is more pronounced than that of the soliton. 

The numerically-determined energies of these solitons have been 
plotted in figure [2j as have the energies of straight solitons determined 
from (57). It can be seen that the energies of Skyrme-Faddeev soli- 
tons determined using full numerical simulations are lower than those 
determined using (57). For larger periods, the difference is small and 
can be attributed to numerical error (the energies determined from 
(57) are more accurate than those determined using full numerical 
simulations). For smaller periods, the difference is more pronounced, 
and occurs because the numerically-determined solutions are slightly 
buckled. 

The elastic rod model gets the broad shape of the energy curve 
right, and energies are predicted with fairly good accuracy for periods 
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in the range [5, 10]. The energy match for smaller periods is not so 
good, however. This is not surprising, because the Skyrme-Faddeev 
model has many more degrees of freedom than the elastic rod model. 
When an elastic rod is compressed beyond P « 5, its energy increases, 
because the winding density 7 becomes large. On the other hand, 
when a soliton is compressed beyond P ~ 5, its energy stays roughly 
constant. The soliton is able to maintain a low energy by changing 
its profile function f(r): figure [4] shows in particular that the soliton 
becomes very narrow at small periods. Elastic rods do not have a 
degree of freedom analogous to the profile function, and this explains 
why rods approximate solitons poorly at low periods. We conclude 
that elastic rods model Q = 1 solitons well qualitatively, but the 
quantitative match is good only for a range of periods. 

5.3 Charge 2 




Figure 5: Energies of solitons and rods with Q = 2 

In figure [5j we have plotted the energies of elastic rods with Hopf 
degree Q = 2 as a function of the period P. For large enough periods, 
the only state in the elastic rod model is a straight rod (29). As 
the period decreases the energy decreases, and at the critical period 
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Figure 6: Elastic rods with Q = 2 and P = 15, 12, 10, 8, 6, 4. 



P = 12.0 the straight rod buckles to form a helix with M = 1. The 
energy of the helix continues to decrease with the period until P = 2.49 
is reached, at which point the non-intersection constraint influences 
the shape of the rod and its energy begins to increase. Table (40) 
indicates that there exists a helix with M = 2, but this always has 
a higher energy than the M = 1 helix and has not been plotted. A 
sample of energy-minimising elastic rods have been depicted in figure 

EES 

Figure [7] displays pictures of energy-minimising Skyrme-Faddeev 
solitons for a range of periods. There is an excellent match with the 
elastic rods displayed in figure 5.3 , including a buckling transition from 
a straight soliton to a helix at a period P ~ 12. 

The energies of the Skyrme-Faddeev solitons have been plotted in 
figure [5j along with the energy of a straight soliton determined from 
( |57[ ). There is a small discrepancy between the energies determined 
from ( |57[ ) and those determined from full numerical simulations, and 
this can be attributed to numerical error. For small periods the ener- 
gies of the numerically-determined solitons are significantly lower than 
those of the straight solitons, reflecting the fact that the solitons in 
the numerical simulations are buckled. 
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Figure 7: Solitons with Q = 2 and P = 15, 12, 10, 8, 6, 4. 



It is clear from figure [5] that the energies of Skyrme-Faddeev soli- 
tons and elastic rods with Q = 2 are in remarkably good agreement. 
We conclude that, for Q = 2, elastic rods are a good model of Skyrme- 
Faddeev solitons, both qualitatively and quantitatively. 



5.4 Charge 3 



In figure |5.3| we have plotted the energies of elastic rods with Hopf 
degree Q = 3 as a function of the period P. For large enough peri- 



ods, the only state in the elastic rod model is a straight rod (29). As 
the period decreases the energy decreases, and at the critical period 
P = 18.03 the straight rod buckles to form a helix with M = 1. The 
energy of the helix continues to decrease with the period, and at the 
critical period P = 13.46 the helix buckles again to form a kinked con- 



figuration, as described in subsection 4.4 This kinked configuration 



cannot be constructed analytically, and has instead been constructed 



using the numerical methods described in subsection 4.5 The ener- 
gies of the numerically-obtained rods dip below the helix energy when 
P < 13.46, indicating that buckling has occurred. We have also plot- 
ted the M = 2 helical state, whose energy is slightly greater than that 



of the M = 1 helix and the kinked configuration. From (40) we see 
that there exist in addition helices with M = —1,3 but their energies 
are much greater than the other states and have not been plotted. 
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Figure 8: Energies of elastic rods with Q = 3 as a function of P 

The elastic rod model predicts, therefore, that there should exist 
kinked Skyrme-Faddeev solitons with Q = 3. In order to test this 
prediction, we have simulated Skyrme-Faddeev solitons with Q = 3 
and P = 6. In fact, we found both a kinked soliton and a helical 
soliton. Both of these are depicted in figure [9j along with the kinked 
elastic rod. The kinked soliton has energy E$f = 2.285 and the 
helix has energy E$f = 2.282. Unfortunately the energies of the two 
solutions are within numerical accuracy, so we cannot conclude which 
has the lowest energy. The helix was obtained by starting with a 
helical initial condition. In order to obtain the kinked configuration 
we started with a configuration obtained from the kinked elastic rod, 
using a construction presented in pp. 

6 Conclusions 

We have investigated minimum-energy configurations in the Skyrme- 
Faddeev and elastic rod models on R 2 x S 1 with Hopf degrees 1, 2 and 
3. We have found a good agreement between the two models, both 
qualitatively and often quantitatively. For all charges and sufficiently 
large periods, the minimum-energy configuration in both models is a 
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Figure 9: Solitons and rods with Q = 3 and P = 6. The elastic rod is on 
the left, the kinked soliton solution is in the middle and the helical soliton 
solution is on the right. 



straight rod (or soliton). Buckling occurs as the period is reduced. 

For Hopf degree 1 , the straight rod (or soliton) buckles slightly to 
form a helix in both models. For Hopf degree 2, a much more visible 
buckling occurs and again the minimum energy configuration in both 
models at small periods is a helix. Elastic rods with Hopf degree 3 
undergo two successive buckling transitions, passing through a helix to 
form a kinked configuration at low periods. The kinked configuration 
also appeared in simulations of the Skyrme-Faddeev model, although 
more numerical work is needed to determine whether or not it has a 
lower energy than a helix. 

Our results show that the elastic rod model is a reliable description 
of Skyrme-Faddeev solitons. They also demonstrate the utility of this 
model: without it, we would not have found the kinked configuration 
at Hopf degree 3. Solitons in the Skyrme-Faddeev model are notori- 
ously difficult to find numerically, particularly at high Hopf degree, 
and it is hoped that the elastic rod model will provide a useful tool 
for tackling this problem. 

More generally, our results demonstrate a clear link between two 
models from field theory and elasticity theory. They motivate the 
search for elastic phenomena in other field theories. 
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